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Abstract. A scattering approach to entanglement in mcsoscopic conductors 
with independent fermionic quasiparticles is discussed. We focus on conductors 
in the tunneling limit, where a redefinition of the quasiparticle vacuum transforms 
the wavefunction from a manybody product state of noninteracting particles to a 
state describing entangled two-particle excitations out of the new vacuum. The 
approach is illustrated with two examples (i) a normal-superconducting system, 
where the transformation is made between Bogoliubov-de Gennes quasiparticles 
and Cooper pairs 1 , and (ii) a normal system, where the transformation is made 
between electron quasiparticles and electron- hole pairs |21 15|. This is compared 
to a scheme where an effective two-particle state is derived from the manybody 
scattering state by a reduced density matrix approach. 
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1. Introduction. 

Over the last decade, entanglement has come to be viewed as a possible resource for 
various quantum information and computation purposes. The prospect of scalability 
and integrability of solid state quantum circuits with conventional electronics has led 
to great interest in the investigation of entanglement in solid state systems. A broad 
spectrum of proposals for generation, manipulation and detection of entanglement in 
solid state systems is given in this volume. Of particular interest is the entanglement 
of individual quasiparticles in mesoscopic conductors. Phase coherence is preserved 
on long time scales and over long distances, allowing for coherent manipulation 
and transportation of entangled quasiparticles. Moreover, individual quasiparticles 
are the elementary entanglable units in solid state conductors and investigation of 
quasiparticle entanglement can provide important insight in fundamental quantum 
mechanical properties of the particles and their interactions. 

Recently, a number of proposals for creation of entanglement in mesoscopic 
systems based on scattering of quasiparticles have been put forth El El E] El E3 
El El El ID] Q2| • 0m main interest here is to discuss a central aspect of such 
systems operating in the tunneling regime, namely the role of the redefinition of the 
vacuum in creating an entangled two-particle state. To this aim we consider the 
two original proposals, Refs. £Q and [21 EJ, were the ground state reformulation was 
discussed. The emphasis in these works was on entanglement of the orbital degrees of 
freedom P^, the discussion however applies equally well to spin entanglement. In 
the first proposal ,T|, we investigated a normal mesosocopic conductor contacted 
to a superconductor. The superconductor was treated in the standard mean-field 
description, giving rise to a Bogoliubov-de Gennes scattering picture with independent 
electron and hole quasiparticles. It was shown that Andreev reflection at the normal- 
superconducting interface together with the redefinition of the vacuum can gives rise 
to an entangled two-electron state emitted from the superconductor into the normal 
conductor. Second, Beenakker et al and later the authors investigated a normal 
conductor in the quantum Hall regime. It was shown that the scattering of individual 
electron quasiparticles together with the redefinition of the vacuum can give rise to 
emission of entangled electron-hole pairs from the scattering region. 

In the present paper, first a general framework for entanglement of independent 
fermionic quasiparticles in mesoscopic conductors is presented. The role of the system 
geometry and the accessible measurements in dividing the conductor into subsystems 
as well as defining the physically relevant entanglement is emphasized. We then 
consider a simple, concrete example with a multiterminal beamsplitter geometry 
and derive the emitted manybody scattering state. Two different approaches to the 
experimentally accessible two-particle entanglement are discussed. First, a general 
scheme for arbitrary scattering amplitudes based on a reduced density matrix approach 
is outlined and then applied to the state emitted by the beamsplitter geometry. The 
orbital entanglement of the reduced state is discussed in some limiting cases. Second, 
for the system in the tunneling limit, we show how an entangled two-particle state is 
created by a redefinition of the vacuum. The redefinition leads to a transition from a 
single-particle to a two-particle picture, with the wavefunction being transformed from 
a many-body state of independent quasiparticles to a state describing an entangled 
two-particle state created out of the redefined vacuum. Based on this discussion, a 
detailed investigation of the entanglement in the systems in Refs. P E] is presented. 
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2. Entanglement in mesoscopic conductors. 

The concept of entanglement appeared in physics in the mid nineteen thirties [13| as 
a curious feature of quantum mechanics giving rise to strong non-local correlations 
between spatially separated particles. The non-local properties of entanglement 
contradicted a commonly held "local, realistic" view of nature and led Einstein, 
Podolsky and Rosen (EPR) to conclude, in their famous paper ^3], that quantum 
mechanics was an incomplete theory. With the inequalities of Bell |15j . presented 
three decades later, it became possible to experimentally test the predicted non- 
local properties of entangled pairs of particles. Since then a large number of 
experiments have been carried out, predominantly with pairs of entangled photons 
[161 1171 I18j . where a clear violation of a Bell Inequality has been demonstrated, 
providing convincing evidence against the local realistic view of nature. To date, 
however, no violation of a Bell Inequality with electrons has been demonstrated. 

During the last decade, the main interest has turned to entanglement in the 
context of quantum information processing . It has become clear that entanglement 
can be considered as a resource for various quantum information tasks, such as 
quantum cryptography |20|. quantum teleportation |21| and quantum dense coding 
22 . The notion of entanglement as a resource naturally led to the question of 
how to quantify the entanglement of a quantum state. A considerable number 
of measures of entanglement have been proposed to date [231 I24j . ranging from 
describing abstract mathematical properties of the state to quantifying how useful 
the state is for a given quantum information task. It is mainly the prospect of 



Figure 1. Picture of orbital (upper) and spin (lower) entanglement scheme. From 
the source S, an orbitally entangled, e.g. \&o) = i /v / 2(|l)A|l)_B + |2}a|2)b), or 
spin entangled, e.g. s) = l/V^d T)a| I) b — \ 1)a\ T)b) state is emitted, with 
one particle propagating towards A and one towards B. At A and B, the particles 
are modulated, e.g. via single qubit rotations parameterized by angles </>a and 
cj>B- The particles are then detected in electronic reservoirs + and — . 



quantum information processing in solid state conductors that has motivated the 
recent interest in entanglement in mesoscopic conductors. In this work we do however 
not try to answer the ambitious question what quantum information tasks can be 
performed in mesoscopic systems. We do also not consider any particular measure of 
entanglement, instead we focus the discussion on the entangled quantum state. Given 
the state, the entanglement, quantified by ones measure of liking, can in principle 
be calculated. Therefore we are interested in the basic first step in a quantum 
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information processing scheme, namely the creation, manipulation and detection of 
quasiparticle entanglement. A large number of implementations of this first step have 
been proposed, typically a certain mechanism is suggested that leads to emission 
of entangled particles from a "source" S. The particles propagate out to spatially 
separated regions A and B, where they are manipulated and detected. Two schematics 
of generic systems for orbital and spin entanglement are shown in Fig. ^ 

Very recently, in a number of works [HElElIHEIEldlHIEl HDl CU entanglement 
in systems of quasiparticles has been investigated within the framework of scattering 
theory. These proposals are of particular interest because working with independent 
particles allows for a complete characterization of the emitted many-body state for 
arbitrary scattering amplitudes. Moreover, the notion of entanglement without direct 
interaction between the quasiparticles, a known concept in the theory of entanglement 
25 , has appeared puzzling to members of the mesoscopic community. The perception 
in the mesoscopic physics community has started to change only with the apperance of 
Ref. [2]. This makes it important to thoroughly analyze the origin of the entanglement. 
Here we contribute to such an analysis, focusing on the creation of entangled two- 
particle states due to redefinition of the vacuum. For comparison, an approach 
based on the reduced two-particle density matrix is discussed as well. Although not 
investigated here, it is probable that the two types of states turn out to be of different 
"usefulness" for quantum information processing, making a detailed comparision of 
interest. We start by stating some important general properties of entanglement and 
comment on their application to entanglement in mesoscopic conductors. 

i) The entanglement of a state in a given system depends on how the system 
is formally parted in to subsystems [20] ■ If the system is considered as one entity, 
i.e. with all quasiparticles living in the same Hilbert space, one has to consider the 
question of entanglement of indistinguishable particles fZT\ EH]. F° r a mesoscopic 
system of independent fermionic quasiparticles, the ground state is given by a product 
state in occupation number formalism, i.e. the wavefunction is a single Slater 
determinant. Apart from the correlations due to fermionic statistics the particles show 
no correlations and the state is not entangled. However, if the system is considered 
as consisting of several spatially separated subsystems, one can pose questions about 
the entanglement between spatially separated, distinguishable quasiparticles living in 
the different subsystems. The latter is typically the case in the proposed mesoscopic 
systems (see Fig. H), which are naturally parted into the source S and the two regions 
A and B. The entanglement between two particles, one in A and one in B, is in many 
situations nonzero. Quite generally, the physically relevant partitions into subsystems 
are defined by the system geometry and the possible measurements one can perform 
on the system |251l3b] . 

ii) The entanglement that can be investigated and quantified is determined by 
the possible measurements one can perform on the system. Although a wide variety 
of measurements in mesoscopic conductors in principle can be imagined, the most 
commonly measured quantities are currents and current correlators, i.e. current noise 
[^Ei]- In particular, in several recent works it was proposed to detect entanglement 
via measurements of current correlators [H El 01 El El El El El CHI IIH 113 ESI E31 This 
leads us to focus the discussion on entanglement detectable with current correlation 
measurements. All current cross correlation measurements, needed to investigate 
spatially separated particles, are to date measurements of second order correlators. 
Importantly, the second order current correlators are two-particle observables and can 
thus only provide direct information about two-particle entanglement. We thus limit 
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our investigations here to the two-particle properties of the state |35| . 

Considering the typical mesoscopic entanglement setup in Fig. EJ the points i) and 
ii) lead us to discuss the entanglement between two spatially separated particles, one in 
A and one in B, detectable via correlations of currents flowing out into the reservoirs. 
Such a bipartite system is also the most commonly considered one in investigations 
of entanglement. Importantly, the entanglement can be both in the orbital pQ as well 
as in the spin degrees of freedom, the relevant entanglement depends on the proposed 
system geometry which can be designed to investigate orbital El El 13 or spin 
@J El El EH EH EH EH1 E3 entanglement. We also note that our measurement based 
scheme excludes all types of occupation number or Fock-space entanglement, e.g. the 
linear superposition of two particles at A and two at B [which in an occupation number 
notation can be written as the Fock-space entangled state 10)^12)^ + 12)^10) a]- There 
is thus no need to enforce additional constraints |36l IHT) on the state to exclude the 
physically irrelevant Fock-space entanglement. 



3. System and scattering state. 



To clearly illustrate the basic principles, we present here the formalism for independent 
quasiparticles in the most elementary mesoscopic system possible. We focus the 
discussion on the orbital entanglement which has the advantage that it can be 
manipulated and detected E] El with existing experimental techniques. However, 
for completeness, spin information is retained throughout the discussion. 




Figure 2. Schematic of the entanglement setup. The source region S (dashed 
red box) consists of two beamsplitters connected to four reservoirs 1, 2, 3 and 
4. The regions A and B (dashed blue boxes) consist of one beamsplitter and 
two reservoirs each. The orbital modes |1) and [2) out of the source region are 
displayed. Reservoirs 1 and 2 are kept at a potential V, the remaining reservoirs 
are grounded. 

The system is shown in Fig. El In the source region, two single mode reflectionless 
beamsplitters are connected to four electronic reservoirs 1,2,3 and 4. A voltage 
bias V is applied to reservoirs 1 and 2, while reservoirs 3 and 4 are kept at zero 
bias. Quasiparticles injected from reservoirs 1 and 2 scatter at the beamsplitters and 
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propagate out towards regions A and B. In A and J5, the quasiparticles are scattered 
at another pair of beamsplitters (i.e. a local single qubit rotation) and then detected 
in four reservoirs A+, A—, B+ and B— , all kept at zero potential. We emphasize that 
all components of this system can be realized in a conductor in the quantum Hall 
regime 

In this work we are interested in the entanglement of the state emerging from 
the source region, before the quasiparticles reach the beamsplitters in regions A and 
B. The manipulation and detection processes taking place in A and B are thus not 
investigated. The two orbital modes |1) and |2), with labels 1 and 2 denoting from 
which reservoir the particles are emerging (see Fig. [2J , constitute the orbital two-level, 
or pseudo-spin, system. Working within the quasiparticle scattering approach [35], we 
introduce operators al n(T (E) creating electron quasiparticles at energy E with spin a, 
incident from reservoirs m = 1, 2, 3 and 4. The energy is counted from the Fermi level 
of the unbiased reservoirs. The transport state of the system at zero temperature is 
given by 

l*>= II 4^E)4 L (E)a^(E)4 L (E)\0) (1) 

0<E<eV 

where |0) is the quasiparticle vacuum, a filled Fermi sea at energies E < 0. This 
describes filled, noiseless streams of electrons emitted from reservoirs 1 and 2, 
propagating towards the first pair of beamsplitters. To describe the state after 
scattering at the beamsplitters, we introduce operators bAm'a-'{E) and bBm'cr'{E) for 
particles originating from reservoir m', propagating from the beamsplitters towards 
regions A and B respectively. The operators bAm'a'(E) and bBm'a'(E) are related to 
a>mcr(E) via the scattering matrices of the beamsplitters. For simplicity we consider 
identical beamsplitters, taken independent on E and <j, giving (suppressing spin and 
energy index) 

b A i\ = (r t' \ ( ax \ ( b A2 \ = f r t' \ ( a 2 
b Bi V t r' { a 3 r { b B2 I t r' J \ a A 



(2) 



We note that only static scatterers are considered here, the discussion could however 
straightforwardly be extended to time dependent scatterers, recently investigated in 
the context of entanglement in Refs. |38l I39| . 

To simplify the notation, a collective quantum number n = {1 |> 1 it 2 T> 2 4} 
denoting both orbital mode and spin is introduced. The state describing particles 
propagating out towards A and B can then be written 



l*) = 

0<E<eV,n 



n Kn(^)+^k(^)]io) (3) 



Due to the scattering, the outflowing streams are noisy and the properties of the 
particles flowing towards A and B can be investigated via noise measurements. Here 
we are interested in the entanglement between two particles in the outflowing streams, 
one towards A and one towards B. Importantly, the state |^) in Eq. J3J is a many- 
body state, it describes a linear superposition of different number of particles at A 
and B, ranging from zero to in principle infinity. Since only two-particle entanglement 
is considered, one thus needs to deduce the two-particle properties of the state |^). 
Below we present two different approaches to do this, giving rise to two different 
quantum states with in general different entanglement. 
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4. Reduced two-particle density matrix approach. 

We first discuss a general approach, applicable for arbitrary scattering amplitudes. 
Starting with the properties of the correlators of outflowing currents towards A and 
B, we note that the cross correlators (in the most general situation) are determined 
by averages of the type [321 

(b\ n (E)b Am (E')b Bk (E")b Bl (E'")) 

= {b\ n {E)b Bk {E")b Bl {E"')b Am {E')) cx p% m (E, E>, E", E">) (4) 

using in the second step the anticommutation relations for fermionic operators. 
The term p nm ki(E, E' , E" , E'") is by definition an element of the reduced two- 
particle density matrix (in the energy basis). The current correlators are thus fully 
characterized by the reduced two-particle density matrix 

p= ]T J dEdE'dE"dE"'p k n l m (E, E', E", E'") 

nmkl 

x bl n (E)b Bk (E'')\0)(0\b Bl (E''')b Am (E'). (5) 

Accordingly, the entanglement potentially detectable via noise measurements is the 
entanglement of the reduced two-particle density matrix. Clearly, this is the situation 
for any two-particle observable. It should be emphasized that the two-particle density 
matrix physically describes the correlations of two particles out of the streams flowing 
towards A and B, leaving all other particles unobserved. This is qualitatively different 
from a projection of a two-particle state out of the full many-body state, where only the 
components of the state containing exactly two particles are selected. It is shown below 
that these two different ways of extracting a two-particle state out of a many particle 
state can give rise to different states, and consequently to different entanglement. 

Even with the object of interest confined to the reduced two-particle density 
matrix, a full characterization of the entanglement is cumbersome since p is in general a 
mixed state and the two particles live in infinite dimensional Hilbert spaces spanned by 
{E, n}. In some simple situations (see e.g. below), the density matrix in Eq. © can be 
written as a direct product of the density matrices for the different degrees of freedom, 
e.g. p = po ® ps ® Pe, where the subscripts 0,S and E denote orbital, spin and 
energy respectively. This allows one to independently characterize the entanglement 
with respect to the different degrees of freedom. In particular, the orbital subspace in 
Fig. [SJand in e.g. Refs. ^EUS] as well as generically the spin 1/2 space are two-level 
systems or qubits, giving rise to a system of two coupled qubits, well studied in the 
entanglement literature (see e.g. Ref. |4"U]1 

However, in the general case, with arbitrary scattering amplitudes, the reduced 
two-particle density matrix can not be written as a direct product, i.e. p 
po® ps® Pe- What can then be said about the entanglement? Of particular interest 
is the entanglement detectable via zero frequency current correlators, generally the 
quantity investigated in experiments. The zero frequency limit effectively projects the 
operators in Eq. |@J to the same energy . For scattering amplitudes independent 
on energy on the scale of eV, the energy argument in Eq. (0) can be suppressed, 
giving a reduced density matrix 

P = Pn m b An b Bk \0)(0\b Bl b Am , p% m cx {b\ n b Bk b B ib Am ). (6) 

nm,kl 
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The density matrix p contains information about the orbital and spin parts of the state 
only. Interestingly, following the opposite approach and considering the quasiparticle 
(E > 0) correlator for coincident times, closely related to the electronic counterpart 
j3] to the joint detection probability introduced by Glauber |41) in quantum optics, 
one finds the same reduced density matrix p as in the zero frequency limit. Similar 
results for short time current correlators have been obtained in Refs. [51 l9l 1101 111) . 

The reduced density matrix p for the system under consideration (see Fig. is 
evaluated from Eqs. QJ, J2J) and © to be 



(lo ® Is - F (8 F s ) , F = 
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(7) 



expressed in the orbital basis {|1)a|1)b, |1)a|2)b, |2)a|1)b, |2)a|2)b} and the spin 
basis {| T)a| T)s, | T)a| |)b, \ 1)a\ T)b, | |)a| 1)b}- Here 1 is the 4 x 4 unit matrix. 
Note that since the particles in A and B are distinguishable there is no need for anti- 
symmetrization and we can use the notation e.g. «AiT a B2.i |0) = |1)a|2)b8>| T)a| |)b- 
Formally, the first term in Eq. J7J) results from the direct pairing of the operators in 
the bracket in Eq. (0 and is given by the direct product of the single particle density 
matrices at A and B. The second term results from the exchange pairing. Since the 
direct term by definition does not describe any nonlocal correlations (it is a separable 
state with respect to A and B), it is the exchange correlations which are responsible for 
the entanglement. For the current correlators, it should be noted that the first term 
determines the product of the averaged currents, while the second term determines 
the irreducible current correlators, i.e. the correlators of the current fluctuations, the 
noise. We also note that p is independent on scattering amplitudes, a consequence of 
the assumption of identical, spin independent beamsplitters. 

As is clear from Eq. J2J), even though we consider a simple geometry with spin- 
independent scattering, the reduced density matrix p is not a direct product between 
orbital and spin part, i.e. p ^ po® Ps- Many observables are however not sensitive to 
the spin degree of freedom, as is the case for cross correlators between total currents 
I = /-)-+ 1 1 , typically the cross correlators investigated in mesoscopic conductors. The 
effective orbital density matrix accessible via current correlators is then obtained by 
tracing p over the spin degree of freedom, giving 

po = tr s [p] = \ (21 - F) = \ (1 + 2|* > (*o|) (8) 
6 6 

with Y$>o) = 1/V2 [|1)a|2)b — |2)a|1)b]- This state, an example of a Werner state 
|42| . can via suitable local transformations be written on a separable form 031 ^] 
with respect to A and B and is consequently not entangled. We note that the same 
holds for the reduced spin density matrix ps, obtained by tracing over orbital degrees 
of freedom. 

A different situation occurs if one considers a spin-polarized system, as was done 
e.g. in [21 El ■ I n this case the density matrix is purely orbital, obtained from Eq. © 
by suppressing the spin notation. For the system in Fig. [21 we obtain 

Po = l(l-F) = \9o){*o\, (9) 

an orbital singlet, i.e. a maximally entangled state, again independent on scattering 
amplitudes. This result can be understood by considering the energy and spin 
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Eq. Q appropriate under the stated assumptions. The corresponding outgoing state 

1 c 



" 2 b A1 b A2 + t 2 b B1 b B2 + rt (b A1 b B2 - b\ 2 b BX 



(10) 



The first two terms describe two particles at A or two at B, while the last term 
describes one particles at A and one at B. As is clear from Eq. ©, only the last term, 
which is just Y&o), contributes to po- Importantly, for the two-particle state in Eq. 
<|10[1 . the reduced density matrix approach gives the same result as projecting out the 
part of the state which contains one particle at A and one at B. Both procedures are 
thus equivalent to a post-selection of entanglement, as originally discussed in quantum 
optics 021 ( see e -g- Ref. for a discussion for fermions). Various issues of projection 
and post-selection were recently discussed in a number of works on entanglement in 
mesoscopic conductors with arbitrary scattering amplitudes [31[7|EI1E] (f° r a related 
discussion, see also [H]). 

It is interesting to note the clear difference between an orbital state obtained by 
tracing p over the spin degrees of freedom and the orbital state in a spin polarized 
system. The difference can be attributed to the fact that only spins of the same spieces 
are nonlocally correlated, i.e. contribute to the entanglement of the state. Detecting 
(without spin resolution) two particles, one at A and one at B, the probability of 
obtaining two identical spins is only one half, reducing the entanglement of the state 
to zero. 



5. Tunneling limit, vacuum redefinition. 

A qualitatively different approach to the characterization of the emitted state can be 
taken in the limiting case of a tunneling system, as was done in Refs. ^ElEj ■ Consider 
the state in Eq. In the tunneling limit, t <C 1, we can expand the product to first 
order in t as 



II [ rb An(E)+tbl n (E) 

0<E<eV,n 

...[rb An (E)+tb Bn (E) 
[rb An (E')+tb Bn (E' 



b An> ( E ) 



tb Bn ,(E) 



[rb\ n {E")+tb Bn {E") 



0<E<eV,n 



where the last term contains fe^n operators at all energies and n's except {E',n'}, 
in the same order as in the line above (we consider a continuous spectrum). 
Using the property that bAn' (E^b^^E')^) — |0) and that the operator products 
bAn'{E')b\ n ,{E') and bBn'{E')^bA n '{E') commutes with all &-operators at different 
energies E ^ E' or n ^ n' , we can write the state in Eq. ® in a suggestive way by 
inserting bA n '(E')b' An ,(E') in front of |0) and reordering the operators, giving 

l*> = II b An(E)\0) 

0<E<eV,n 
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Figure 3. Pictures for particles flowing out from a beamsplitter in the source 
region, a) Single particle scattering picture. A transmission of a particle in the 
filled stream incident on the beamsplitter leads to creation of an electron flowing 
out towards B and a missing electron, i.e. a hole, in the streams flowing out 
towards A. b) Electron-hole pair creation picture, the single particle scattering 
event creates an electron-hole pair on top of the redefined vacuum, with the 
electron moving out towards B and the hole towards A. 



+ t r dE'J2blAE')b Ar AE') n 

-'0 n^fi 1 ^ 



0<E<eV 



l+t 



eV 



IJ b\ n (E)\0) (12) 

. 0<E<eV,n 

The first term describes filled streams of quasiparticlcs flowing out towards A. The 
second term, with a small amplitude t -C 1, describes the same filled streams with 
one electron missing and in addition a second electron propagating out towards B 
(in mode 1 or 2). It is thus very natural to incorporate the filled, noiseless stream 
of quasiparticles flowing out towards A into a new vacuum |0), the second term then 
describes an electron-hole excitation out of the redefined vacuum. Formally, we can 
write the new vacuum in terms of the old vacuum |0) as 



16) = LI b ln(E)\0) 



(13) 

0<E<eV,n 

The new ground state is thus a filled Fermi sea at energies E < in B and at E < eV in 
A. The fcrmionic operators describing excitations out of the new vacuum are defined 

as 



c\^E') = ±b An (E), c\ 



XE) = bl n (E) 



(14) 



where n denotes a spin opposite to n, +(— ) is for spin up (down) in n and E' = 
eV — E. The transformation in Eq. (|14fl is thus equivalent to an inverse Bogoliubov 
transformation. We can then write the wavefunction in Eq. (|12|) as (reintroducing 
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|0) + |*) 



(15) 



with 



!*>=*/ dE [4 W (^)4 T j(^') " 4 Tj -(^) c A|i(^)] |0> (16) 



■eV 



3 = 1,2 



The wavefunction I'l') describes a two-particle excitation, an electron-hole pair, out 
of the redefined vacuum |0). The redefinition of the vacuum thus gives rise to a 
transformation from a picture with a many-body state of independent particles to a 
picture with two-particle excitations out of a ground state. 

There are a number of important conclusions to be drawn from the result above, 
and to be put in relation to the result for the reduced density matrix approach: 

(i) The state is wavepacket-like, i.e. it consists of a sum of electron-hole pairs 
at different energies, a detailed characterization of a similar wavepacket state emitted 
in a normal-superconductor system is given by us in Ref. 07] • As pointed out in 
Ref. pQ, the average time between two subsequent wavepackets is much longer than 
the "width" of each wavepacket, i.e. subsequent entangled pairs are well separated 
in time, in contrast to the reduced two-particle state in the large transparency limit. 
This temporal separation probably makes the entangled state created by a redefinition 
of the vacuum more useful for quantum information processing, due to the possibility 
of addressing the individual entangled pairs. We remark that the formal procedure 
used above to calculate the state in Eq. l(T^|) is essentially identical to the one for 
calculating the wavefunction of the two-photon state emitted in a parametric down 
conversion process in optics |48|. 

(ii) In a first quantization notation, we can write the wavefunction of the 
excitation 



This state is a direct product of the orbital, spin and energy parts of the state, in 
contrast to the general situation for the reduced two-particle state. Both the orbital 
and the spin state are maximally entangled Bell states. The spin state is a singlet 
as one would expect of an excitation out of a spinless groundstate, created by spin- 
independent scattering. 

(iii) As was emphasized in Ref. j^j, the redefinition of the vacuum is possible 
only in fermionic systems, i.e. it relies on the existence of a filled Fermi sea such that 
a removal of an electron below the Fermi energy creates a hole quasiparticle. This 
is further emphasized by the fact that the new groundstate, just as the initial one is 
noiseless. 

(iv) The correlators between electron currents are simply related to the correlators 
between electron and hole currents as (Ar(t)AI h (t')) = -(AI e (t')Ar(t)). 
Consequently, the electron-hole correlators are experimentally accessible and the 
electron-hole entanglement is, in line with the discussion above, a physically relevant 
object to study. 

Importantly, the redefinition of the vacuum and the transformation to an electron- 
hole picture can be performed for an arbitrary transmission. However, in this case the 
resulting state describes a superposition of different numbers of electron- hole pairs. To 




® [| 1)a\Vb-\Va\1)b]®\E') a \E) b 



(17) 
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obtain a two-particle state, one has to calculate a reduced electron-hole density matrix 
along the same line as in Eq. (J5J, i.e. replacing the electron operator a\ n with the 
quasiparticle operator c^ Afl etc. Performing such a calculation in the low transparency 
limit, one obtains as expected with \ty) given by Eq. (|1 T|) . It is however not 

possible from the reduced density matrix approach to conclude whether the emitted 
state is a true two-particle state or a reduced two-particle state. 

For arbitrary scattering amplitudes, to quantitatively compare the reduced 
density matrix approach for electrons and holes to the results for electrons discussed 
above, we consider the simplest situation with low frequency correlators, i.e. all c- 
operators at equal energy, and a spin polarized conductor. The reduced orbital density 
matrix for the system in Fig. [21 is then given by 

2(1 +T) [ T1 + 2R \Vo)(*o\\, (18) 

with |* ) = 1/V2[\1) a \1)b + |2)a|2) b ] and the scattering probabilities R= 1—T = 
\r\ 2 . Interestingly, in contrast to the density matrix in Eq. 10, the density matrix Pq 
depends on the scattering probabilities. The state Pq" is a Werner state, entangled for 
R> T. In the limit T -C 1, one has pQ = | ^o) {^O I : an orbital Bell state, maximally 
entangled. Moreover, away from the tunneling limit, the two-particle entanglement in 
the electron-hole picture is smaller than the entanglement in the electron picture [the 
state in Eq. @ is maximally entangled]. This is in agreement with the findings of 
Lebedev et al ^0]> wn0 investigated the conditions for a violation of a Bell Inequality 
for a scatterer with arbitrary transparency, comparing the electron and the electron- 
hole approaches. 

The state p e Q can be understood by considering the state in Eq. (|10|l after 
redefining the vacuum and transforming to an electron-hole picture, giving 



eh 

Po 



+ t2c Bl c B2 c Al c A2 + rt ( c Al c Bl + c A2 c B2) 1°) ( 19 ) 



The last term is just \^o) while the second term, describing four particles, i.e. two 
electron- hole pairs, also contributes to the reduced density matrix in Eq. (|18fl . it gives 
rise to the term Tl. Interestingly, performing a projection of 1^) onto a state with only 
two particles, one in A and one in B, one obtains the maximally entangled state \^o), 
since the terms with two electron-hole pairs are discarded. The projection approach 
thus, in the large transparency limit, overestimates the entanglement detectable via 
current correlations. 

To illustrate the relevance of the vacuum redefinition, we now discuss two different 
systems where this was investigated. 



6. Normal-Superconducting entangler. 

We first consider the case of a normal-superconducting system, to large extent 
following Ref . . As in all existing works on entanglement in normal-superconducting 
systems @] 03 El EE E3 EH E3 EU 155] - we consider the superconductor in 
the mean-field description. The mean field Hamiltonian is bilinear in fcrmionic 
operators and is diagonalized by a Bogoliubov transformation, giving rise to a new 
set of independent electron and hole like quasiparticles. Importantly, although 
the microscopic mechanism for superconductivity is interaction between electron 
quasiparticles, in the mean field description it is again possible to find a picture 
with noninteracting quasiparticles. For a noninteracting normal conductor connected 
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to a superconductor, the whole system can thus be treated within a single particle 
scattering approach to the Bogoliubov-de Gennes equation [SSj. Consequently, the 
approach above to the entanglement for systems of independent quasiparticles can be 
applied. 




Figure 4. Schematic of the conductor. Left: Mesoscopic conductor contacted to 
a normal reservoir (blue shaded) and a superconducting reservoir (red shaded). 
Right: Paths of scattering particles at energies well below the superconducting 
gap. An incoming electron from the normal reservoir has the amplitude r ee to be 
reflected as an electron and r he to be reflected as a hole. For incoming holes, the 
corresponding amplitudes are r eh and r . 



To illustrate the basic principle of two-particle emission, we first consider the case 
of a two-mode (to connect to the orbital discussion above) normal-superconductor 
system shown in Fig. At the normal-superconductor interface, for energies well 
below the superconducting gap, scattering occurs either as Andreev reflection or 
as normal reflection. Consider the situation with a negative bias — eV applied at 
the normal reservoir while the superconducting reservoir is kept at zero potential. 
Counting the energy from the superconducting chemical potential us, we consider 
operators a en [E) and a\ n (E), creating electron and hole quasiparticle excitations 
respectively at energy E incident from the normal reservoir. The collective quantum 
number n denotes as above orbital mode and spin. The state of the system at zero 
temperature is given by 

\*)ns= II a L(EM (20) 

0<E<eV,n 

describing a filled stream of holes injected from the normal reservoir, where |0) 
is the quasiparticle vacuum. The operators a en {E) and ah n (E) are related to 
operators b en (E) and bh n (E) for quasiparticles propagating back from the normal- 
superconducting interface towards the normal reservoir via the scattering matrix 
(again taken independent on E and n) as 



b e 

b h 



r eh 
Jih 



(21) 



(22) 



The state |$)jvs can then be written in terms of the 6-operators 
Wns= II [r eh bL(E)+r hh b{ n (E)]\0) 

0<E<eV,n 

Considering the tunneling limit, the Andreev reflection amplitude r eh <C 1 and one can 
proceed in exactly the same way as in Eqs. i|ll[l and (|12fl to arrive at the wavefunction 
to leading order in r eh as 

l + r eh f dE'J2blAE')b hn iE>) [] 
Jo , 



NS 



(EM 



0<E<eV,n 
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(23) 

The first term describes a filled stream of hole quasiparticles flowing back towards the 
normal reservoir. The second term, with a small amplitude r eh <C 1, describes the 
same filled stream with one hole missing and in addition a second electron propagating 
back. We can thus proceed as above and incorporate the filled, noiseless stream of 
hole quasiparticles flowing out towards the normal conductor into a new, redefined 
vacuum |0), given by 

10)= II 6 U£)|0> (24) 

0<E<eV,n 

The fermionic operators describing excitations of the new ground state are given by 
the Bogoliubow transformation 

4 (E) = b\ n {E), cU-E) = ±b hn (E) (25) 

where n denotes the spin-flip in the electron-hole transformation with +(— ) for spin 
up (down) in n. Note that the c-operators are just standard electron operators. A 
missing hole in the filled hole stream is thus just an electron with opposite spin and 
energy compared to the superconducting chemical potential us- We can then write 



a) 



eV 



b) 



Figure 5. Schematic of a normal-superconducting interface, a) Bogoliubov-de 
Gennes picture, a filled stream of holes are incident on S from N. An Andreev 
reflection leads to creation of an electron above the superconducting chemical 
potential fig and a missing hole in the back- flowing hole stream, b) Cooper pair 
tunneling picture, a Cooper pair tunnels from S to N, leading to a pair of electrons 
in N, on top of the ground state. 

the wavefunction in Eq. (|12[) as, reintroducing spin and orbital notation 

\*}ns= |0)jv S + |*) (26) 

with 



|*) = r 



eh 



[ dE^2 kV^hV-^) - c U E ) c U- E )] i°) ( 27 ) 

J ° 3 = 1,2 1 ' J 
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The wavefunction |\&) thus describes a two-particle excitation, an electron pair or 
Cooper pair, out of the redefined vacuum |0). The redefinition of the vacuum thus 
gives rise to a transformation from a picture with a many-body state of independent 
particles to a picture with two-electron excitations out of a ground state. We point 
out that this approach, shown schematically in Fig. provides a formal connection 
between the scattering and the tunnel Hamiltonian approach (see e.g. Refs. |47l 149) 1. 

To connect this two-particle emission to orbital entanglement, we consider as 
a concrete example our proposal in Ref. jj. The system geometry is shown in 
Fig. El a multiterminal normal conductor connected via tunnel barriers to a single 
superconductor and further to four normal reservoirs. The two regions at A and B, 
with an electronic beamsplitter and two reservoirs respectively, constitute the two 
subsystems. For details of the proposals we refer the reader to Ref. and Ref. 
57 . To first order in tunnel barrier transparency a pair of electrons is emitted on 
top of the new ground state, at interface 1 or 2. This leads to a state describing a 




Figure 6. Orbital normal-superconductor entangler: A single superconductor 
(S) is connected to four normal arms via two tunnel barriers 1 and 2 (thick black 
lines). The arms are joined pairwise in beam splitters A and B and end in normal 
reservoirs labeled + and — . After Ref. Q- 



linear superposition of pairs at 1 and at 2. Taking into account that the electrons 
can be emitted either towards A or towards B, we have the emitted state given by 
|*)ivs = |*i)ats + \^2)ns, where 

r eV 

r 

j=l,2r]=A,B 



J _♦_ 1 Q __ A D 



\^2)ns = r e 



"V3l 
eV 



(E)c^(-E) |0) 



dE 



eV 



cV£)4 u (-£) - <4u(£)4ir(-i?) 



+ ^(E)cl n (-E) - c\ 2l (E)cl n (-E) |0), 



(28) 



The state |\E'i) ns describes a superposition of two particles at A and two at B. This 
state does however not contribute to the noise correlators or, as is clearly the case, to 
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the two-particle density matrix describing one particle at A and one at B. The state 
\^z)ns however describes orbitally (as well as spin-entangled) electron "wave packet" 
pairs, with one particle at A and one at B. This is the entangled state detected by the 
noise. 

7. Normal state entangler. 

We then turn to a normal state entangler working in the Quantum Hall regime. Such 
a system was introduced by Beenakker et al in Ref. and then later considered in |7j . 
In both Refs. [2] and [7], two parallel edgestates connected via a nonadiabatic, edge 
channel mixing scattering region was considered. In Ref. [3] we instead considered a 
topologically different Quantum Hall system, a Hanbury Brown Twiss geometry (or 
Corbino geometry) with only single edge-states and quantum point contacts. This 
highly simplifies the experimental realization of the proposal. The system is just a 
proposal for a physical realization of the system shown in Fig. [21 discussed in detail 
above, and we therefore keep the discussion short. 

The system is shown in Fig. for details we refer the reader to Ref. 
A spin-polarized edge state is considered. We take the transmission and reflection 
probabilities at the point contact C to be Tc = 1 — Rc = T and at D to be 
To = 1 — Re, = R. After scattering at C and D, the state I*) consists of two 
contributions in which the two particles fly off one to A and one to B, and of two 
contributions in which the two particles fly both off towards the same detector QPC. 
Consider now the case of strong asymmetry R -C 1, where almost no electrons are 
passing through the source QPC's towards B. Performing the reformulation of the 
ground state and the transformation ta an electron-hole picture, we can directly write 
the full state I*) to leading order in VR as I*) = |0) + VR\^>), with 



Due to the redefinition of the vacuum, we can interpret the resulting state as 
describing a superposition of " wavepacket" -like electron- hole pair excitations out of 
the new vacuum, i.e. an orbitally entangled pair of electron-hole excitations. This is 
just equivalent to the result in Eq. l(T7l) . 

8. Conclusions 

We have presented a scattering approach to entanglement in mesoscopic conductors 
with independent fermionic quasiparticles. The role of the system geometry and 
accessible measurements in defining the relevant entanglement was discussed. As 
a simple example, a multiterminal mesoscopic conductor with spin-independent 
scattering was investigated in detail, deriving an expression for the manybody 
scattering state emitted by the conductor. The focus was on orbital entanglement, 
accessible with present day experimental technics. 

Two different approaches to a two-particle state were considered. The main focus 
was on conductors in the tunneling limit, where entangled two-particle states arise 
from a redefinition of the quasiparticle vacuum. The redefinition of the vacuum leads 
to a transition from a picture with a many-body state of independent particles to a 
picture with two-particle excitations out of a new ground state. This was compared 
to the entanglement of an effective two-particle state, obtained by a reduced density 




■ eV 



(29) 
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2 




3 

Figure 7. Normal-state orbital entangler: a rectangular Hall bar with inner and 
outer edges (thin dashed lines) and four quantum point contacts (grey shaded) 
with transparencies Ta , Tb , Tc and To . Contacts 2 and 3 are sources of electrons 
(a voltage eV is applied against all other contacts which are at ground) . Electrons 
follow edge states (thick black lines) in the direction indicated by the arrows. An 
Aharanov-Bohm flux $ penetrates the center of the sample (shaded). After Ref. 
0. 

matrix approach, applicable for conductors with arbitrary scattering amplitudes. 
Moreover, the qualitative difference between the reduced two-particle density matrix 
approach and a projection scheme was discussed. 

We showed that the two different approaches, redefinition of the vacuum state and 
two-particle state reduction, in general give rise to different two-particle states and 
consequently to different entanglement. This was illustrated by investigating in detail 
the orbital entanglement for the simple mesoscopic conductor, focusing on the spin- 
polarized regime. In the tunneling limit, we applied our approach with the redefinition 
of the vacuum state to the systems proposed in Refs. 1 , a normal-superconducting 
heterostructure, and in [3], a conductor in the Quantum Hall regime. This showed the 
qualitative as well as quantitative similarities between the entangled states emitted in 
the two systems. 
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